Inhomogeneities and impurities in a dense one- dimensional Rydberg lattice gas 
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We consider a dense one-dimensional laser-driven Rydberg lattice gas with perfect nearest- 
neighbor blockade. The ground state of this system can be found analytically in certain parameter 
regimes even when the applied fields are inhomogeneous in space. We will use this unique feature to 
investigate the effect of an impurity - introduced by the local variation of the laser parameters - on 
the correlations of the many-body ground state. Moreover, we explore the role of a staggered laser 
field which alternates from site to site thereby breaking the sublattice symmetry. We demonstrate 
that this technique, which can be applied experimentally, reveals insights into the role of long-range 
interactions on the critical properties of a Rydberg gas. Our work highlight novel possibilities for 
the exploration of many-body physics in Rydberg lattice gases based on locally tuneable laser fields. 



Lattice gases of Rydberg atoms are currently in the 
focus of intense research as a platform to investigate 
strongly-correlated quantum systems. A host of theo- 
retical works has studied and analyzed many-body phe- 
nomena such as the formation of crystalline structures as 
well as their melting under the influence of quantum fluc- 
tuations or dissipation pTfTB]. These investigations have 
revealed detailed insights into the equilibrium properties 
and non-equilibrium phenomena of strongly interacting 
many-body quantum systems. 

On the experimental side there has recently been much 
progress. While initial experiments had revealed first ev- 
idence for the strong and coherent interactions present 
among Rydberg atoms pTH22] a very recent experimen- 
tal breakthrough has enabled the unambiguous identi- 
fication of interaction effects manifesting themselves in 
strong spatial correlations among Rydberg atoms excited 
from a large lattice [23] . This new class of experiments 
permits not only the direct visualization of Rydberg cor- 
relations but also the manipulation of external fields on 
the scale of individual lattice sites. This opens up new op- 
portunities for the study of many-body phenomena and 
the exploration of phase transitions in Rydberg lattice 
gases. 

In this work we perform a theoretical investigation of 
the ground state of a one-dimensional Rydberg lattice gas 
in the presence of inhomogeneities which are induced by 
spatial variations of external laser fields. On the one hand 
these spatial variations are intrinsic to any experiment 
and therefore an understanding of their effect is of tech- 
nical relevance. On the other hand deliberately imposed 
inhomogeneities can be used to experimentally tackle im- 
portant problems in condensed matter physics, such as 
the investigation of impurity phenomena and the detailed 
study of statics and dynamics of many-body quantum 
systems exposed to symmetry breaking fields. We show 
that despite the presence of such spatial inhomogeneities 
an approximate analytical description of the Rydberg lat- 
tice gas in terms of a frustration-free Hamiltonian can 
be obtained. We will use this feature to perform a de- 
tailed investigations of the effect of an impurity and al- 
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FIG. 1. Schematics of the one-dimensional lattice with spac- 
ing, a. Each site contains a single atom whose ground state 
1^)^ is coupled to a Rydberg nS'-state |t);. via a laser param- 
eterized by a site-dependent Rabi-frequency Qk and detuning 
Afc. (a) Rydberg lattice gas with a single impurity present 
on the j-ih site. The impurity atom is irradiated by a laser 
(orange) of different parameters (^^imp, Aimp) compare to the 
rest of the system (blue, and parameterized by Qgys, Asys). 
(b) Alternating lasers are introduced to investigate the effect 
of sublattice symmetry breaking fields. On odd (even) lattice 
sites, a laser parameterized by Qodd (^even) and Aodd (Aeven) 
is used to excite atoms to Rydberg states. 



ternating laser fields which break the sublattice symme- 
try of the Rydberg gas Hamiltonian. This work comple- 
ments current studies on Rydberg gases which are almost 
exclusively conducted on homogeneous systems [5H3 M- 
[TT| 116] . It illustrates the usefulness of the frustration- free 
description and sheds light on the critical properties of 
the ground state of a dense Rydberg lattice gas. 

We consider a system in which atoms are held in a 
deep one-dimensional optical lattice with inter-site spac- 
ing a (see Fig. fl]) and with a single atom per site, i.e., 
the external atomic degrees of freedom form a Mott- 
insulator. To model the internal dynamics we employ the 
well-established two-level pseudo-spin description: The 
electronic ground state of the atom located at the k-th 
lattice site is denoted by |^)^. This state is coupled to 



a Rydberg n5'- state, denoted by |t)/., through a laser of 
Rabi frequency fl^ and detuning A/g, both of which can 
be site-dependent. Atoms in Rydberg nS'-states interact 
via the van-der-Waals interactions V{r) = Cq/t^ where 
Cq is the dispersion coefficient and r the interatomic dis- 
tance. With these prehminaries we can now formulate 
the Hamiltonian of this many-body system on a lattice 
with L sites: 

L L L 

HRyd = Y^ ^kcrl + ^ AkUk + Yl y\k-m\nknm' (1) 

k k m>k 

Here Uk = (l + cr|)/2 = |t)/e(t| is the projector on the Ry- 
dberg state and a^ and cr^ are Pauli spin matrices. The 
interaction strength is parameterized by the constants 
Vk = Ce/{akf. 

In the following we will show that this Hamiltonian 
permits an approximate analytic solution over a certain 
range of parameters. To this end we generalize a con- 
struction put forward in Refs. [7 and consider a regime 
in which the laser Rabi-frequency is much weaker than 
the interaction among nearest neighbors, i.e. Vi ^ Q^- 
Here the simultaneous excitation of neighboring atoms 
to Rydberg states is energetically forbidden. This sup- 
pression is made manifest in the Hamiltonian i^Ryd by 
projecting it onto the subspace in which excited neigh- 
boring spins are absent, H = V i^Ryd V, using the projec- 
tion operator V = Yli^i^ — n^n/e+i). Only the first term 
of i^Ryd is affected by this projection. It transforms ac- 
cording to V CF^V = Pk-i c^k -^k+ii where P^ = 1 — n/e is 
the projector on the electronic ground state of the /c-th 
atom. 

We proceed by neglecting, for the moment, interac- 
tions beyond next-nearest neighbors, i.e. we set V^ ^ 
for /c > 2. Within these approximations the projected 
Hamiltonian of the Rydberg gas assumes the form H = 



- Y.k ^k^k + ^RK, with 
L 

Hrk = Y^^^Y1 ^kPk-1 {<yl + £.kPk + £>k^nk) Pfe+1, 
k k 

(2) 

provided that the laser parameters and the interaction 
strength obey the relationship: 



nk = V2i^\ Ak = nk{Ck 



Hk) 



This can be checked by direct inspection utilizing the fact 
that one is working in the subspace V where nkUk+i = 0. 
Eqs. (|3| describe a region in parameter space where the 
projected Rydberg Hamiltonian is a so-called frustration- 
free or Rokhsar-Kivelson Hamiltonian [24, 25 . Here the 
ground state ll^n})^ which depends on the set of param- 
eters {^n}, is annihilated by all local Hamiltonians hk 
such that hk \{£,n}) =0. It can be constructed analyti- 



cally and its explicit form reads 
1 



|{^n}) = 



V^iM 



exp 



X^aP,_ia+Pfc+i)|0), (4) 

k J 



where Z^^^y is a normalization constant, |0) = \m - ■ - i) 
is the spin vacuum and a^ = (a^ + ^^D/^ ~ lt)/eUI- 
The relationship between the parameters ^^ and the laser 
parameters, i.e. 1^/c, A^ , is obtained by eliminating V2 
from Eqs. (|3|: 



a = 



6 




(5) 



The square of ^k can be interpreted as local fugacity since 
the probability for a spin on the k-th site to be in the 
state |t)/j. is proportional to ^|. This can be directly 
seen by expanding Eq. Q in the spin product basis. 
In the following we will therefore refer to ^k as fugacity 
parameter. 

Let us now make use of the analytical knowledge of the 
ground state and investigate a Rydberg lattice gas in the 
presence of an impurity. In this scenario, which is de- 
picted in Fig. [ija), all atoms but the j-th are irradiated 
by the same laser field characterized by l^sys and Agys 
and the corresponding fugacity parameter ^gys^ satisfy- 
ing the relation ([3|. The impurity is introduced on site 
j by applying here a laser field with in general different 
parameters given by Ai^p and l^imp- The corresponding 
fugacity parameter is then ^imp, which in an experimen- 
tal setting can be achieved by single site addressing as 
experimentally demonstrated in Ref. [23] . 

In Fig. [2] we show the spatial density distribution cal- 
culated for a lattice with periodic boundary conditions 
and length L = 41. The impurity is introduced at site 



J 



21 and its fugacity parameter is set to ^[ 



imp 



5. 



The fugacity parameter of the system atoms ^gys is var- 
ied on the vertical axis. For small ^gys a peak of high 
density is visible at the site of the impurity which is ex- 
pected since ^imp > ^sys and thus, according to Eq. Q, 
the excitation probability at site j, is enhanced. When 
^sys is increased the overall density increases and density 
oscillations in the vicinity of the impurity indicate the 
onset of stronger and stronger correlations in the ground 
state. Indeed one can show that the correlation length of 
(3) the system spins is given by Ic = ^sys [7 • This length is 
plotted as the cyan diagonal lines in the figure showing 
that the impurity indeed determines the state of the sys- 
tem atoms located within a distance Ic. In the vicinity 
of 6mp = ^sys the density modulations vanish as here the 
system becomes translationally invariant. 

When the fugacity parameter ^sys is increased towards 
infinity - a limit which is achieved when (^gys approaches 
zero - the correlation length diverges. Here the ground 
state Q of the homogeneous system, i.e. if the impu- 




FIG. 2. Rydberg density for a lattice with 41 sites and an 
impurity placed at site j = 21. The fugacity parameter of the 
impurity is given by ^imp = 5 (black dash) and the fugacity 
parameter of the remaining spins ^sys is varied from to 10. 
With increasing ^sys the system's correlation length Ic (plotted 
at the cyan diagonal lines) increases. The many-body state 
of the system atoms becomes more strongly correlated and 
the presence of the impurity affects the state of more and 
more distant atoms. Note that for ^sys = 6mp the density is 
homogeneous in space due to the translational symmetry of 
the system. 



rity was absent, becomes a superposition of the two anti- 
ferromagnetic states I tit It •••) and |ttttt •••)• Moreover, 
one can show that the excitation gap of the Hamiltonian 
([2| closes and in fact the impurity breaks the sublattice 
symmetry of the Hamiltonian Q by forcing the ground 
state into one of the anti-ferromagnetic states with a 
down- spin at the impurity site. 

To investigate the symmetry breaking region in pa- 
rameter space in detail we introduce a staggered laser 
field. In an experimental setting this would be done by 
choosing laser parameters such that they give rise to a 
fugacity parameter ^odd (^even) at all odd (even) lattice 
sites, c.f. Fig. [ijb). For the upcoming analysis it is 
convenient to define the difference (_d = ^odd — ^even and 
the sum ^^ = ^odd + ^even of the two sublattice fugacity 
parameters. To investigate the effect of the symmetry 
breaking field let us now specifically study the expecta- 
tion value of the Rydberg density on the odd sublattice 
^odd = X]/c=odd ^k- Usiug the grouud state Q one finds. 



(^odd) 



U-^7' 



V(i+^^"')(i+a, 



(6) 



For small ^~^ this function predicts a transition between 
two states in which Rydberg atoms predominantly oc- 
cupy the odd/even sublattice which takes place when 
the difference between the sublattice fugacity parameters 
vanishes ^^ = [see Fig. Isja)]. This is expected since 
for ^1 = and ^rf = both |ttttt •.•) and |;tttt •.•) 




FIG. 3. (a) Mean density of the odd sublattice as a function 
of ^d- Here, with ^s = 5, we have plotted the analytical result 
given in Eq. Q in black and the numerical result obtained 
from diagonalizing Hamiltonian (IT]) in red (with circles) with 
L — 14. The latter shows a significantly steeper switching 
of the sublattice populations at ^^ = 0. (b) Susceptibility 
Xodd(^s = 5, ^d) for different lattice sizes: L = 10 (red circles), 
L = 12 (blue triangles), and L = 14 (black squares). The 
data suggests a divergence of the susceptibility at ^d = in 
the limit of large lattice sizes L. 



are ground states and any non-zero value of £^d will fa- 
vor one over the other. Note, that according to Eq. (IgI) 
this switching between the sublattices is not sharp as the 
susceptibility, i.e. the slope Xodd(6,^d) = S(^odd)/9^d, 
saturates at a value 1/2 at the "transition point" {^^ = 
0,^~^ = 0} [see Eq.([6|]. Since a tiny perturbation to the 
system would lead to a symmetry breaking, e.g. an impu- 
rity, one would naturally expect the transition occurring 
at ^(i = to be sharp rather than a crossover. Therefore, 
although the frustration-free Hamiltonian ^ excellently 
describes the Rydberg gas along the curve parameter- 
ized by Eqs. ([3| as shown in Ref. [7 , it is very ques- 
tionable whether this Hamiltonian and Eq. ^ faithfully 
describes the actual sublattice occupation of the ground 
state of the Rydberg gas Hamiltonian ^ at the "tran- 
sition point" . The suspicion is confirmed by numerically 
calculating (riodd) in the ground state of i^Ryd- This 
data is shown as the one with red circles in Fig. IsFa) and 
clearly displays a significantly sharper transition. More- 
over, as shown in Fig. [SJ^b), one can anticipate a diverging 
behavior of the susceptibility Xodd{£,s^£,d) with increasing 
lattice sizes. 

This strongly suggests that {^d = O^^^^ = 0} is a 
critical point of the Rydberg gas Hamiltonian (fTj) which 
is not captured by the ground state of the frustration- 
free approximation ^. To investigate the nature of this 
point we perform a scaling analysis by using the results 
shown in Fig. |3|b). Expressing the susceptibility as 
Xodd ^ \^d\~^ with critical exponent 7, in the vicinity 
of ^d = O5 we extract the critical exponent by fitting 
log I Xodd I as function of log|^rf| linearly within a linear 
scaling region [26 and subsequently determining the gra- 
dient. For ^5 = 10 and L = 18, within an appropriately 
chosen scaling region, the critical exponent is found to 
be 7 ~ 1.76 ± 0.05. This result suggests that the second 
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FIG. 4. Stylized phase diagram of the homogeneous Ryd- 
berg lattice gas Hamiltonian (IT]) with Vt^ — Vt and A^ = A. 
The curve parameterized by Eqs. (|3| approaches in the limit 
^ ^ CO the critical point {Q = 0, A = — 3V2} which is lo- 
cated between a classical crystalline phase with Rydberg den- 
sity 1/2 and one with density 1/3. Inset: In our analysis we 
approach the critical point {^d — 0,^7^ = 0, V3 = 0} from 
different directions: One experiences a 2D-Ising type transi- 
tion by approaching the critical point from finite ^d, while one 
experiences a AFM-PM transition by approaching from finite 

V3. 



order transition belongs to the 2D-Ising universality class 

At the critical point we have ^^ = and therefore 
^odd = ^even which means that the fugacity parameter 
is site independent, (,k = £,- To understand the limits of 
the validity of the frustration- free Hamiltonain ^ and to 
gain further insights into the critical point it is therefore 
instructive to have a closer look at the phase diagram of 
the Hamiltonain (IT]) in the homogeneous case. A stylized 
version of this diagram is depicted in Fig. |4| Here one 
sees that the curve, which is parameterized by Eqs. (|3| 
hits the point {Q = 0, Ac = —3^2}, i-e. ^ ^ 00, where 
the phase boundaries of a classical phase with Rydberg 
density 1/3 and one with density 1/2 coalesce. An inspec- 
tion of the frustration-free Hamiltonian ([2| shows that in 
the limit ^ ^ 00 this approximate Hamiltonian does not 
energetically discriminate between configurations of the 
form Itltltl •••) and Itlltii •••) due to the lack of the 
next-next-nearest neighbor interaction Vs. This explains 
why its ground state (p| fails to describe properties, such 
as the sublattice density of the Rydberg lattice gas, in 
the vicinity of the critical point. 

Having in mind that the presence of Vs can affect the 
critical behavior, let us finally investigate in more de- 
tail the role of long-range interactions. To this end we 
amend the frustration-free Hamiltonian ([2| by adding 
next-next-nearest neighbor interactions: Hs = i^RK + 
^3 S/c^/c^fc+3- While in the previous inhomogeneous 
case, the critical point was approached from finite ^d and 
growing ^ = ^s/2, we approach it now on the V3 axis. 
This is depicted as the inset of Fig. [4| We continue 



by deriving an effective Hamiltonian emerging from Hs 
in the limit ^ ^ 00 following the procedure outlined in 
Refs. |9ll27]: We introduce as new degrees of freedom the 
spin blocks 1 =ti and (D =tii- States composed out of 
a concatenation of these spin blocks are strictly degen- 
erate under the frustration-free Hamiltonian ([2| when 
^ ^ 00. This degeneracy is broken by the next-next- 
nearest neighbor interaction V3, and, at finite but large 
^, due to the presence of virtual transitions of the type 
tii^iti- I^ order to derive an effective Hamiltonian in 
this parameter regime one identifies the above-mentioned 
spin block 1 as fictitious hard core particle with creation 
operator 6j., and the block as corresponding fictitious 
hole. Following Refs. [Ql, l28] and introducing the spin 
operators S^ = (hk + 6^)/2, S^ = i(hk - hl)/2 and 
S^ = h\hk — 1/2 one finds that the effective theory for 
large ^ is given by the Heisenberg spin- 1/2 XXZ-model 
in a magnetic field. 



H^ 



= E 



Wx'^x 



+1 



Q3 CJ + 1 



) + fisi-"' 



/i 



Si 



with 3 = -1/3 and /i = 2/3 + ^ V3. With respect to 
the pseudo-spin degree of freedom, i.e. |t)^, and ||)^, 
this model exhibits a critical transition between an anti- 
ferromagnetic (AFM) and a paramagnetic (PM) phase 
at /i = (5 + 1 [29], i.e. when ^ V3 = 0. Hamiltonian Hs 
has thus a critical point at {V3 = 0,^~^ = 0}. This 
corroborates the previous analysis of the numerical data 
displayed in Fig. |3] which suggested a second order tran- 
sition of the Rydberg lattice gas Hamiltonian (fTl) to oc- 
cur at {^d = O?^^^ = 0}- III fact due to the finite 
nearest-neighbor and next-next-nearest neighbor interac- 
tion which is inevitably present in practice the Rydberg 
gas will never be strictly at the critical point. However, 
due to the smallness of V3, stemming from the r "^-scaling 
of the van der Waals interaction, observable exhibits a 
scaling behavior (as for example shown in Fig. ^ which 
will be interesting to be explored experimentally pursu- 
ing the route in Refs. [2ll3Q]. 

In conclusion, we have presented an analysis of the 
static properties of a dense inhomogeneous Rydberg lat- 
tice gas based on a frustration- free Hamiltonian. This ap- 
proach allowed us to analyze in detail the effect of an im- 
purity on the correlation properties of the system. More 
importantly, the introduction of a symmetry breaking 
field allowed us to identify a critical point and to study 
the influence of long-range interactions on the critical 
properties of the system. In the future it will be in- 
teresting to extend this study to other inhomogeneous 
situations: For example, the generality of Hamiltonian 
([1]) permits the exploration of disordered systems. More- 
over, the fact that near the critical point the physics is 
effectively described by the XXZ-model might enable the 
experimental implementation and study of impurities im- 
mersed in Luttinger liquids along the lines of Ref. [31] . 
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